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Abstract 

The Weyl-Wigner-Moyal formalism is developed for spin by means of a correspondence between 
spherical harmonics and spherical harmonic tensor operators. The analogue of the Moyal expansion 
is developed for the Weyl symbol of the product of two operators in terms of the symbols for the 
individual operators, and it is shown that in the classical limit, the Weyl symbol for a commutator 
equals i times the Poisson bracket of the corresponding Weyl symbols. It is also found that, to the 
same order, there is no correction in the symbol for the anticommutator. 
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For a particle with linear momentum p and coordinate q, the Weyl-Wigner-Moyal formal- 
ism jl-3] leads to a fascinating correspondence between classical and quantum mechanics, 
as it provides a mapping between quantum mechanical operators and classical dynamical 
variables defined as functions on phase space. When the operator is the density matrix, 
the phase space function is known as the Wigner function. Because this function fails to be 
nonnegative in all circumstances, it cannot be taken to be a true probability distribution. So 
while it does nothing to solve the knotty problem of interpreting quantum mechanics, this 
formalism is the closest one can get to turning quantum mechanics into a purely classical 
statistical theory, and it also finds practical application, especially in quantum optics. 

It is desirable to have a similar formalism for spin degrees of freedom both for intrinsic 
reasons, and for developing the semiclassical limit. In their low energy states, many atoms, 
molecules, and molecular ions behave as particles with a fixed magnitude of spin that is 
often large, in which case a semiclassical approach is natural. We have found, for example, 
that the Weyl representation is advantageous in constructing the spin-coherent-state path 
integral [4]. The problem of developing such a formalism has been approached by several 
workers over the years fl-EJ with varying emphases and methods. The problem is solved in a 
formal sense and its abstract mathematical aspects are quite highly developed. At the same 
time, unlike the case of linear position and momentum degrees of freedom, there are few 
results in closed form or of practical value. To our knowledge, it has not even been shown 
that the leading classical behavior of the Moyal bracket as the spin magnitude j — > oo is the 
same as the Poisson bracket. This result is of course to be expected on physical grounds, 
but that is not the same as actually proving it, and the proof is an important physical check 
on the correctness of the formal structure. One can further expect higher order corrections 
in powers of 1/j to aid in developing quantum corrections. While our eventual goal is to 
address this practical problem, in this paper we only present general results. Our approach 
is more succint in some ways than the formally oriented ones, relies more on analysis than 
algebra, and may appeal to physicists for these reasons. 

That a Weyl representation for spin operators has been known for several decades may 
surprise some readers, as it surprised us. We have found modern day papers which state the 
contrary, which along with our own ignorance illustrates the surprisingly common occurrence 
in modern day physics of the same facts being well known to one subcommunity of physicists 
and totally unknown to others. Therefore, we cannot hope too much that our paper will 



lessen this disconnect. 

Our basic working tool is the set of spherical harmonic tensor operators 3^ m (J). We find 
that the P, Q, and Weyl representations (indicated by the superscripts P, Q, and W) of 
yt m (J) are all proportional to Y £m (h). That is, for a particle of spin j, 

( n ) = a je ^m(n), (1) 

where the coefficients are given by 
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Since any operator can be uniquely and prescriptively expressed as a linear combination of 
the 3^m(J)'s, these results constitute an algorithmic solution to the problem of finding the 
P, Q, and Weyl representations for an arbitrary operator. 

Before proving these assertions, we review some basic definitions and known results. 
The objects of interest are operators for a particle of spin j, which are functions of J, the 
vector spin operator, with components J x , J y , and J z that obey the usual commutation rules 
(setting h to 1) [J a , Jp] = ie a ^ 7 J-y, and J ■ J = J% + Jy + J z = j(j + 1). For an arbitrary 
such operator A(J), the P and Q symbols are defined by 

$J(n) = (n|A(J)|n), (5) 
A(3) = ^±1 f rfn|n)^(n)(n|. (6) 
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Here, |n) is the (normalized with (n|n) = 1) spin state with maximum spin projection along 
the direction n. That is, Jn |n) = j\h). In particular, with \ j, m) denoting the simultaneous 
eigenstate of J • J and J z , |z) = \j,j). In Eq. (JHJ), the integral is over all n. It also pays 
to employ stereographic coordinates on the sphere. If 9 and tp denote the usual spherical 
polar coordinates of n, and we define z = tan \6e llf ', then it is a standard result that the 
(unnormalized) state \z) defined by 

\z) = e zJ ~\j,j) = V ( . 23 ) z>~ m \j,m) (7) 
z — ' V 7 — ml 



is on the same ray in Hilbert space as |n). That is, \z) = (1 + \z\ 2 ) 3 ; |n), where the normal- 
ization follows from Eq. (J7J). When A = !V^ m (J), we write simply ^f^ ,W for the symbols. 

We also recall that Eq. (jB} does not define $^ uniquely. The reason is that the matrix 
elements of |n) (n| in the \ j, m) basis are linear combinations of Yg m (n) with £ < 2j. Because 
of the orthogonality of the Yg^s, we may add to any linear combination of Y^ m 's with 
£ > 2j without affecting the integral over n. In this paper we shall fix <3>^ by demanding 
that the coefficients of Yg_ m with £ > 2j vanish. 

Any physically meaningful phase space representation must be linear, real for Hermitean 
operators, equal to the number 1 for the identity operator, and covariant under rotations. 
That is, if the symbol for an operator A(3) is $^(11), and if a rotation takes J into J' and 
n into n', then the symbol for A(J') must be $^(n'). These properties are obeyed by the 
P and Q symbols (modulo the nonuniqueness of the former). The Weyl symbol must also 
obey the key extra demand of traciality j^]. That is, if A and B are two operators, we must 
have 

' tr (AB) = ±. /dn$J(nX(n). (8) 
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The left side of this equation may be seen as a quantum mechanical average, and the right 
side as a classical one over phase space. We denote these as ( ) qm and ( )„. This condition 
determines $^(n) unambiguously. Further, if B, say, is taken as the density matrix, then 
$g (n) is the Wigner function for the system. 

Next, we must define the operators 3^m(J)- We do this in parallel with the spherical 
harmonics Y^ m (n) via the Herglotz generating function 10] for the latter: 
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Here, 
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^(* + iy) + ^(*-iy), (11) 



with A and ( being real. The essential property of a is that a • a = 0. Further, r = |r|, 
and n = r/r. Equation (Q implies the normalization that the integral of |Y£ m | 2 over all n is 
unity. 



Equations © and (jTOl) show immediately that various relationships among the Y im 's have 
direct analogs for the 3^ m 's. We have, for example, 

K/Jn) = (-l) m F,,_ m (n), y\ m {3) = (-l) m ^,- m (J). (12) 

Further, it is apparent that under rotations Y^ m (n) and 3^ m (J) transform identically. 

We now note that 3^ m (J) is the operator analogue not of the surface harmonic Y^ m (n), 
but of the solid harmonic r l Yi m (n) which has the extra factor r l . Hence, we expect all 
three symbols for 3^m(J) to asymptote to j^Y^ m (n) as j — > oo. Second, for spin j, any 
operator is equivalent to a matrix of order 2j + 1. A Hermitean matrix of this order has 
(2j + l) 2 independent parameters, so there can be at most (2j + l) 2 independent 3-Ws. This 
number is exhausted by taking £ up to 2j, and the rotational properties then guarantee that 

y tm (j) = o$t>2j. 

We now find the various $| m 's. For <3>2n we take the expectation value of Eq. (jUJ) in the 
state |n) to obtain 
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We now evaluate the left hand side using the stereographically labeled states {\z)}. Thus, 

(n|e ?a - J |fi) = (1 + \z\ 2 )- 2l {z\e^ 3 \z) 

= (l + \z\ 2 r 2 i{j,j\e** J +e< a - J e zJ -\j,j). (14) 

The Lie algebra of the J a 's allows us to write e z * J + e ( a - J e zJ - — e u-J- e pJz e u+J+ ^ w h ere u± anc [ 
(3 are functions of (, z*, z, and a. We thus obtain (,z|e^ |z) = e™, and (n|e^ a ' J |n) = e 7 - 7 , 
where e 1 = e /3 /(l + \z\ 2 ) 2 . 

To find 7, we exploit the faithfulness of the j = 1/2 representation of SU(2). In this case, 
(a- J) 2 = aiaj(5ij + 2ieijkJk)/4: = 0, so, e^ a J = l + £a-J. Now, (n|J|n) — jh, so for j = 1/2, 
(n|e^ a ' J |n) = 1 + |£a ' n = e7,/2 - Hence, for general j, Eq. (fT3l) may be written as 

(n|e^ J |n> = (l + l( a .h) 2j (15) 
2 i /o„\ H 
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where in the last line we have used the ( e term in Eq. (|9]). Comparing with Eq. (fT3|) . and 
simplifying, we obtain 

which is the same as Eq. ([3]). As expected, ^^(n) i s proportional to 5^ m (n), and it vanishes 
if £ > 2j. Since any operator is determined by just its diagonal matrix elements in the spin 
coherent states, we have rededuced that 3^m = for I > 2j. Alternatively, knowing this 
result, we could have anticipated Eq. ( jTTl) from the connection between solid and surface 
harmonics which ensures that $^(n) ~ j l Y im {h). 

Knowing $^(n), §f m is simple to find. Taking the expectation value of Eq. (jfJJ) in the 
state |n), and noting that |(n|n')| 2 = [(1 + n ■ n')/2] 2 - ? , we obtain 

^-^/^(IlfJ^VKtf). (18) 

It is straightforward to show that 
where 

^' = (2<+1) W -oS'L + i)i - (20) 

and is the Legendre polynomial of order £ We take A = 3^m(J), and feed and $^ 
from Eqs. ([I])-© along with Eq. (fT9l) into Eq. (fT8l) . and invoke the addition theorem and 
orthonormality of the Y^ m 's to obtain 

Reading off aj t from Eq. (fT7j) . we immediately obtain Eq. (j2J) for a^. 

We now find $^(n). We evaluate the trace in Eq. (JSJ) by writing Eq. (JBJ) for A, and 
multiply by B to obtain 

<$^(n«(n)) fi = (A5) qm = <^(n)^(n)) fi . (22) 

Naturally, we can switch the roles of P and Q on the right. If we take 3^m for A and 3Vm' 

<^ m (J)^L'(J)) qm = ^ W) V^', (23) 



for B, we find Eq. (j4]) for when £ = m = m'. More generally, 



which is an operator orthogonality relation. In particular, = + l)] 1 ^ 2 , so that 
for J itself, the Weyl symbol is \fj{j + 1) n. This is a direct deduction of exactly what is 
indicated for the "classical" vector corresponding to J by a large number of indirect quantum 
mechanical arguments. Further, for fixed £, we find that 

*2>(J)« (l + ! + ---)-7%»( fi )> tf->oo). (24) 

The above formulas allow us to obtain a miscellany of other results. First, we can invert 
Eq. (llSp to go from the Q symbol to the P symbol: 

<^(n) = i- j dn'G(n,n')$S(A'), (25) 

where 

c^^e ^^^ P^y (26) 

The cooefficients in this formula are essentially 1/Kj£, so it is somewhat formal, and it would 
be interesting to obtain a result in closed form. 

Second, we can find asymptotic relations between the different symbols in the classical 
limit, by which we mean that we let j — > oo keeping fixed the form of the operator A as a 
function of J. To get <3>^ from e.g., we take Eqs. (J2J) and (jBJ), and expand the ratio in 
powers of 1/j to get 
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Now, $f m = a^Ja%, and £(£ + l)^ m (n) = C 2 Y em (h), where £ = -i(n x V fi ) is the 
angular momentum operator (on phase space, and noi the quantum mechanical Hilbert 
space). When written in terms of £ 2 , the ratio a^Ja% has no £ dependence, and so constitutes 
an operator which applied to will yield for any A. Up to order j -3 , we find 11] 

= (l + + ^^(^ 2 - 2) + - 2 )(^ 2 - 3 ) + • • (28) 

Similarly, 

= (l + ^ + ^ £2 ( £2 " 4 ) + ^ £ ^ £4 - 8^ 2 + 24) + • • • ) $^(n). (29) 

We can also find the Stratonovich-Weyl kernel A J (n) 0,0], which is is an operator valued 
function on the unit sphere with the properties that 

$J(n) = {AAi(&)) A = <<^(n)A^(n)) fi . (30) 



It is straightforward to show that 
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W 

•it- 



(31) 
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instead, A J is given as a matrix in the \ j, m) basis. It is easily seen that A J must 



obey the rules, 



(A'(A)>* = 1, <A'(ft)A'(A')> = 4^(n,n'). 



(32) 



Here, P'(n, n') is the identity kernel on the space of harmonic functions on the unit sphere 
that are of order 2j or less. That is, 



J dh'P(h,h')f(h') = f(h) 



(33) 



for any function /(n) whose spherical harmonic expansion contains no terms with £ > 2j. 
The first property of A J follows almost trivially from Eq. (j3ip . and the second from Eq. (1231) . 

We turn, finally, to determining the asymptotic (j — > oo) behavior of the Moyal product, 
i.e., the Weyl symbol for a product of operators AB in terms of the symbols for the individual 
operators A and B. Using the A- 7 mapping, one finds that js] 



M 3 (h u h 2 ,h 3 ) = (A^nOA^n^A^na)}^, 



(34) 
(35) 



As j — > oo, the kernel Mj is very sharply peaked when rii = n 2 = £13 and nearly zero 
otherwise. Progressively higher moments of Mj are of higher order in 1/j, suggesting a 
Taylor expansion of $^ and $^ . Now, 



(36) 



where Q\ 2 is the angle between n 2 and n 1; and £±,12 is the component of £ orthogonal to 
both hi and ri 2 . We have 

#12 



12^±,12 
6*12 
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(fii x n 2 ) • £ 
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(37) 



sin 9x2 

We thus generate a series in powers of gradients of and Omitting the arguments 
of Mj , we have 

(M i /(n 4 )> = /(n 1 ), (z = 2,3), (38) 



for any /(n). Hence, beyond the leading term, all terms in the Moyal product will entail 
gradients of both and The first non trivial term (of order arises from the 

moment 

(M^cJizfi) = - ( w , 2 ni a n 1/3 + — w e a ^n ly + c5 a/3 , (39) 
\ a ji ) ^ a ji 

where c = (aj^ — j(j + l))/3(a]^) 2 = 0(j -2 ). The first term contributes nought, and the 
third only at order 1/j 2 . The central term leads to 

-U • x £<) = -Lfl • (V fi $T x V fi <), (40) 

2j c 2j c 



where j c = + 1) and we have written n for fii . 

We can find moments of higher order tensors such as r^a 7 ^/??^ in the same way. Odd 
order tensors do not lead to any corrections, and only even orders need to be considered. 
However, even finding all fourth order moments and the corresponding part of the Moyal 
product of order j~ 2 is laborious and quite daunting. So, in summary, what we have found 
up to now is that 

SJk(A) = <(n«(A) + ^fi • (V fi $T x V fi <) + 0(r 2 ). (41) 

The term of order 1/j may be written more physically in terns of the Weyl symbol for the 
commutator [A, B]. If we define j c h = j c as the classical angular momentum vector, we may 
cast it as 

But the right hand is precisely the Poisson bracket of two functions of j c , given that the 
fundamental Poisson brackets are {j ca , Jc/3}pb = Ca/3-yjcy Hence we have shown that to 
leading order as j — > 00, 

*1^«W.^}pb- (43) 
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